Ultracold Bosonic Atoms in Disordered Optical Superlattices 
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The influence of disorder on ultracold atomic Bose gases in quasiperiodic optical lattices is dis- 
cussed in the framework of the one-dimensional Bose-Hubbard model. It is shown that simple 
periodic modulations of the well depths generate a rich phase diagram consisting of superfluid, Mott 
insulator, Bose-glass and Anderson localized phases. The detailed evolution of mean occupation 
numbers and number fluctuations as function of modulation amplitude and interaction strength 
is discussed. Finally, the signatures of the different phases, especially of the Bose-glass phase, in 
matter-wave interference experiments are investigated. 



Quantum phase transitions of interacting bosonic 
many-body systems in disordered lattice potentials have 
been a topic of intense theoretical investigations in the 
past years. A rich variety of possible zero-temperature 
phases has been predicted such as the Mott-insulator 
phase, a Anderson localized phase, and a Bose glass phase 
Gl ||, ||, |J] . Recently, impressive experiments on the tran- 
sition from a superfluid to a Mott-insulator || demon- 
strated that ultracold bosonic atoms in optical lattices 
offer unique possibilities to explore the phase diagram 
of these systems. The degree of experimental control is 
remarkable: The geometry of the lattice potential can 
be designed and specified precisely, even disorder can be 
introduced in a controlled manner || The strength 
of the two-body interaction can be chosen by means of 
Feshbach resonances. Ultimately, the structure of the 
ground state can be examined in detail. This makes them 
a promising candidate to study the competition between 
disorder and interaction experimentally in an unprece- 
dented way. 

In this Letter we discuss the ground state phase dia- 
gram of an interacting bosonic many-body system in an 
optical superlattice, i.e., a quasiperiodic lattice composed 
of disordered unit cells. We study the dependence of the 
ground state on the interaction strength and the amount 
of disorder within the Bose-Hubbard model. 

Consider a gas of N bosonic atoms in an one- 
dimensional lattice potential at zero temperature. For a 
sufficiently strong lattice it is convenient to describe the 
state of the system in a basis of localized "tight-binding" 
wave functions, i.e., the Wannier functions that result 
from a band-structure calculation. We assume that only 
the localized ground state of each lattice well contributes 
and excited vibrational states can be neglected. For a 
system with I lattice sites the many-body state can thus 
be represented in terms of number states |m, ■■■,ni) with 
occupation numbers rij for the individual sites i = 1, /. 
A complete basis of the model space is formed by the 
set of number states \nf,...,nf) (a — 1,...,D) with all 
compositions of the occupation numbers. The dimension 
D = (N + I— l)!/[iV! (7-1)!] of the number basis grows 
dramatically; for I = 8 lattice sites and N — 8 particles 
the dimension is 6435, for / = N — 10 it is already 92378. 



Within the model space any state can be expanded in this 
number basis 



D 



(i) 



with expansion coefficients C a . We introduce creation 
operators and annihilation operators a\ , which create 
and annihilate, resp., a particle in the lowest vibrational 
state at site i, and the corresponding occupation number 
operators = afaj. 

The Hamiltonian of the interacting many-body sys- 
tem in second quantization is the so-called Bose-Hubbard 
Hamiltonian @, |, |, f§ 
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H = [ - J (a\a l+1 + a\ +1 ai) 

»=i 



(2) 
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The first term describes the coupling between neighbor- 
ing sites with a strength J. We use cyclic boundary con- 
ditions, i.e. hopping between the first and the last site of 
the lattice is included. The last two terms of (||) give the 
on-site single-particle energy and the on-site two-body 
interaction with a strength V. Formally, the parame- 
ters are given by matrix elements of components of the 
coordinate space Hamiltonian calculated with the Wan- 
nier wave functions associated with the individual sites. 
The hopping strength J, e.g., is the off-diagonal matrix 
element of the kinetic energy operator calculated in the 
Wannier basis. For an irregular lattice potential the on- 
site energies depend explicitly on the site index i; for 
simplicity we neglect the site-dependence of J and V. 
Interactions between particles at different sites and long 
range hopping are also neglected. 

To determine the ground state of the system we solve 
the eigenvalue problem of the Bose-Hubbard Hamilto- 
nian numerically. The Hamilton matrix in the number 
basis is easy to calculate; the on-site energy and the two- 
body interaction form the diagonal and the hopping term 
generates a few off-diagonal matrix elements. Since the 
Hamilton matrix is very sparse and we are interested in 
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the lowest eigenstates only, an iterative Lanczos-type al- 
gorithm is most efficient to solve the eigenvalue problem. 
This enables us to treat systems with up to / = 12 and 
N = 12 on a standard PC without further approxima- 
tions. Larger systems can be treated by Monte Carlo 
techniques |8[ |9|. Two relevant ground state observ- 
ables we consider in the following are the mean occu- 
pation number fii = {ip\ fii \tp) and number fluctuations 
of = (0| nf \tf)) — (tp\ rii \ip) 2 at the individual sites. 

The interplay between the three terms of the 
Bose-Hubbard Hamiltonian (^) generates a rich zero- 
temperature phase-diagram. Its basic structure can be 
understood by analyzing the contributions of the differ- 
ent terms to the energy expectation value for particular 
states. The off-diagonal hopping term gives a negative 
contribution to the total energy if the state is a super- 
position of many number states. If only the hopping is 
present the coefficients of the ground state (|l|) are related 
to the multinomial coefficients, thus all number states 
contribute and number fluctuations <7j are large. The 
particles can tunnel freely through the lattice and the 
system resembles a superfhiid. 

If a repulsive two-body interaction is included then 
number states with large occupation numbers at indi- 
vidual sites have high energy expectation values. The 
repulsion favors homogeneous distributions of the par- 
ticles over all sites. The competition between hopping 
term and two-body interaction governs one kind of quan- 
tum phase transition present in these systems. For small 
interaction strengths V / J the hopping term dominates 
and the ground state is a superposition of many number 
states — a superfhiid. With increasing V/J those num- 
ber states with large occupation numbers at some sites 
are gradually suppressed because of their large interac- 
tion energy. 

For integer (commensurate) fillings N/I there is a 
unique number state with m = N/I at each site which 
minimizes the expectation value of the two-body inter- 
action. For sufficiently strong interactions the ground 
state is given by just this number state. Therefore num- 
ber fluctuations <7j and also the expectation value of the 
hopping term vanish. This phase in which tunneling of 
the particles is inhibited by the repulsive interaction is 
called Mott insulator phase. Quantum Monte Carlo cal- 
culations for the infinite one-dimensional Bose-Hubbard 
model with N/I — 1 show that the superfhiid to Mott- 
insulator transition occurs at (V/J) C ijt ~ 4.65 [9fl, which 
is confirmed by a renormalization group study [ lCy and a 
strong-coupling expansion (J] . Mean-field models predict 
a much larger critical interaction strength Q [Tl| . 

How does disorder affect the phase diagram of the 
Bose-Hubbard model? Most of the theoretical investiga- 
tions on disorder induced effects concentrate on infinite 
lattices with completely random on-site energies. More 
relevant for cold bosonic atoms in optical lattices are 
quasiperiodic structures, i.e., lattices composed of disor- 
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FIG. 1: On-site energies 
ti for a unit cell of the 
sinusoidal superlattice. 



dered unit cells. As simplest example we discuss an one- 
dimensional optical superlattice with a sinusoidal mod- 
ulation of the on-site energies e^. A unit cell consists of 
1 = 8 lattice sites and the q vary in the interval [— A,0] 
as shown in Fig. [j]. Experimentally this modulation can 
be realized using a superposition of two optical standing 
waves with appropriate wavelengths — a so-called two- 
color lattice. This allows the independent control of the 
tunneling strength J and the disorder amplitude A. The 
sinusoidal modulation already shows the relevant funda- 
mental features but is only the most simple realization of 
disorder, more complex superlattices were already gener- 
ated experimentally ||. 

Let us start from a noninteracting superfluid in a regu- 
lar lattice (A/ J = 0, Vj J = 0) and increases the disorder 
amplitude A gradually. The mean occupation number fii 
of the site with lowest q in each unit cell will increase. For 
sufficiently strong disorder the ground state will be com- 
posed of number states which have non-zero occupation 
only at the deepest well of each unit cell. This mech- 
anism is similar to Anderson localization in an infinite 
system with completely random on-site energies ^] and 
we will adopt the name for simplicity. However, in con- 
trast to the Anderson localized phase in a random lattice 
the ground state is still a superposition of several number 
states with considerable number fluctuations <7j. There- 
fore, this state appears almost like a superfluid state in 
a regular lattice with / times the lattice spacing. 

The disorder induced localization is strongly affected 
by repulsive interactions which drive the system towards 
even distributions of particles over all sites. To study the 
detailed interplay between interactions and disorder we 
diagonalize the Hamiltonian numerically for an isolated 
1 = 8 unit cell. This assumes that the exchange of par- 
ticles between different unit cells can be neglected. By 
direct comparison with two cell calculations this turns 
out to be a very good approximation for V/ J > 0. 

Figure || shows the mean occupation numbers and 
the number fluctuations Ui for the sinusoidal superlattice 
with A/ J = 60 and commensurate filling (/ = 8, N = 8) 
as function of the interaction strength V/J. For vanish- 
ing interaction strength all particles are localized at the 
lowest energy site within the unit cell. However, with 
increasing V / J the mean occupation at this site is re- 
duced rapidly and the particles are redistributed to sites 
with higher on-site energies. This rearrangement hap- 
pens stepwise, i.e., first the two sites with second lowest 



3 




FIG. 2: Mean occupation numbers hi (a) and number fluc- 
tuations Oi (b) for the sinusoidal superlattice with 1 — 8, 
N = 8, and A/ J = 60 as function of V/J. The different rows 
correspond to the individual sites of a unit cell. 

€i are populated. Then at V/J « 10 the population of 
the next two sites increases, and so on. Between suc- 
cessive rearrangements there are extended plateaus (e.g. 
for 20 < V/J < 40) with constant and approximately 
integer mean occupation numbers hi. This characteristic 
sequence of rearrangements and stable plateaus is also re- 
flected in the number fluctuations at shown in Fig. |^(b). 
Within the rearrangement regions number fluctuations of 
those sites which change their mean occupation number 
are large. Within the plateaus all <Ji are small and the 
ground state is almost a pure number state. This region 
of successive rearrangements is called Bose-glass phase 

§!• 

Eventually, at V ~ A a final rearrangement happens 
and all sites have equal hi — 1 and very small number 
fluctuations. This is the transition from the Bose glass 
to the Mott insulator phase where the ground state is a 
pure number state with rii = 1 at each site. 

Figure ^ summarizes the phase-diagram for the sinu- 
soidal superlattice. The contour plot shows the square of 
the largest coefficient C^ ax = max(C^) in the expansion 
(|I|). The dark shadings indicate that the ground state 
is a pure number state. The Mott insulator phase (Mf) 
appears only for V > A. For interaction strengths below 




FIG. 3: Contour plot of the square of the largest coefficient 
Cmax for the sinusoidal superlattice with / = 8, N — 8 as 
function of V/J and A/ J. Dark shadings indicate that the 
ground state is a pure number state (insulators). The labels 
identify the different phases: superfluid (SF), Mott insulator 
(MI), Bose glass (BG), and Anderson localized (AL). 



this boundary two stable configurations within the Bose- 
glass phase (BG) follow. At very small V/J and large 
disorder amplitudes the Anderson localized phase (AL) 
appears. The remaining regions resemble a (disordered) 
superfluid (SF). 

The structure of the Bose glass phase, i.e. the se- 
quence of rearrangements, depends strongly on the par- 
ticular pattern of on-site energies. More complex mod- 
ulations of the ti generate a larger variety of configura- 
tions in the Bose glass phase. The same holds true for 
non-commensurate fillings, fn the sinusoidal lattice with 
N = 9 extended regions emerge, where two sites exhibit 
mean occupation hi = 0.5 or hi — 1.5 in association with 
large number fluctuations. 

How can these structures be observed experimentally? 
An experimentally quite simple approach is the imaging 
of matter-wave interference patterns after release from 
the trap and ballistic expansion Js], [lj| . The crucial quan- 
tity that determines the presence or absence of an inter- 
ference pattern is the phase coherence between atoms at 
different lattice sites. Small phase fluctuations at the 
individual sites, which imply phase coherence, are con- 
nected to large number fluctuations o^. An estimate for 
the phase fluctuations g&a can be obtained from the re- 
lation (70, i — 1/(2(7;) To simulate the density in- 
terference pattern for a given ground state we approxi- 
mate the wave function of an atom localized at a lattice 
site by a Gaussian wave-packet |H| . The superposition 
of the wave packets after free expansion (neglecting in- 
teractions) results in a density interference pattern. To 
account for the phase fluctuations we calculate the inco- 
herent average over typically 10000 sets of on-site phases, 
which are chosen randomly with a Gaussian distribution 
of width cr^^. 

Figure ^ depicts characteristic interference patterns for 
different points of the phase diagram. The superlattice 
character is explicitly included by replicating the unit 
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tive fringe pattern. Eventually, if we enter the Mott insu- 
lator phase only the incoherent bump remains as shown 
in panel (f). To obtain complementary information on 
the spatial density distribution of these disordered states 
Bragg diffraction || |l4|, [l5| would be a useful tool. 

In summary, we have shown that ultracold atomic 
gases in optical superlattices are an ideal system to study 
the complicated interplay between interaction and dis- 
order. As function of the interaction strength one can 
observe the detailed evolution of the ground state from 
an Anderson localized state through a Bose-glass to the 
Mott-insulator, which is accompanied by a characteristic 
vanishing and reappearance of interference fringes. The 
remarkable degree of experimental access to all relevant 
parameters allows a comprehensive study of these quan- 
tum phase transitions. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft, the UK EPSRC, and the EU under the 
TMR network ERB-FMRX-CT-0002. 



FIG. 4: Matter- wave interference patterns for different com- 
binations of V/J and A/ J. The gray curves show the full 
interference pattern of the superlattice; the black curves re- 
sult from folding with a Gaussian profile to mimic a restricted 
experimental resolution. 



cell. This gives rise to the short wave length oscilla- 
tions which reflect the periodicity of the superlattice. The 
black curves are obtained by folding with a Gaussian pro- 
file to mimic a restricted experimental resolution. The 
remaining smooth pattern is largely determined by the 
structure within the unit cells. 

In the absence of interactions and disorder we obtain 
the prominent interference pattern of the superfluid j^] 
with pronounced peaks shown in Fig. |](a). This struc- 
ture dissolves if we enter the Anderson localized phase 
shown in panel (b) for A/ J = 60, depending on the 
experimental resolution one may detect the interference 
pattern of a lattice with / times the fundamental lat- 
tice spacing. If we include interactions V/J > and 
enter the Bose glass phase we observe a characteristic 
vanishing and reappearance of interference fringes which 
is correlated to the number fluctuations shown in Fig. 
^(b). For weak interactions V/J — 10 the three-peak 
structure of the superfluid reappears with broader peaks 
and increased incoherent background due to the reduced 
number fluctuations. In the stable regions of the Bose 
glass phase, e.g. at V/J — 30, number fluctuations of 
all sites are small and the fringes vanish again; only the 
incoherent background remains as shown in panel (d) . In 
the vicinity of rearrangements certain lattice sites regain 
large number fluctuations and cause the reappearance of 
interference fringes. Figure ^(e) depicts an example for 
V/J = 42, where every second lattice site has consider- 
able fluctuations (compare Fig. ||) which create a distinc- 
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